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Abstract 

Tropical differential equations are introduced and an algorithm is 
designed which tests solvability of a system of tropical linear differential 
equations within the complexity polynomial in the size of the system 
and in its coefficients. Moreover, we show that there exists a minimal 
solution, and the algorithm constructs it (in case of solvability). This 
extends a similar complexity bound established for tropical linear sys¬ 
tems. In case of tropical linear differential systems in one variable a 
polynomial complexity algorithm for testing its solvability is designed. 

We prove also that the problem of solvability of a system of tropical 
non-linear differential equations in one variable is AP-hard, and this 
problem for arbitrary number of variables belongs to NP. Similar to 
tropical algebraic equations, a tropical differential equation expresses 
the (necessary) condition on the dominant term in the issue of solvability 
of a differential equation in power series. 

Keywords: tropical differential equations, polynomial complexity solving 


Introduction 

Tropical algebra deals with the tropical semi-rings Z + of non-negative integers 
or Z + U { 00 } endowed with the operations {min, +}, or with the tropical semi¬ 
fields Z or Z U { 00 } endowed with the operations {min, +, —} (see e. g. [6j, 

0 . i). 

A tropical linear differential equation is a tropical linear polynomial of the 
form 

min{ap' ) + x^\ a} (1) 

hJ 
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where the coefficients a, aP e Z + U{oo}, and a variable xp is treated as ” j-th 
derivative of x t := xP". 

For a subset Si C Z + we define the valuation 

Val Si {{j > 0}) := Val Si ({xP}j> 0 ) : Z+ ->■ Z+ U {cx)} 

of variable Xi as follows. For each j > 0 take the minimal s G Si (provided 
that it does exist) such that s > j and put ValsfJ) s — j: in case when 
such s does not exist put Valsjj) oo. We use a shorthand 

Val Sl ,..., Sn := Val Sl X ■ ■ • x Val Sn : Tf (Z+ U {(»})". 

Observe that if W; is a power series in t with the support {f s , s G 5)} then 
Valsi(j) is the order ord t {XP) at zero of the j-th derivative X^\ 

We say that Sj,..., S n is a solution of the tropical linear differential equa¬ 
tion (QJj if the minimum min^-jap' 1 + Vals^j ), a} is attained at least twice or 
is infinite (as it is accustomed in tropical mathematics i. 0). The latter is 
a necessary condition of solvability in power series in t of a linear differential 
equation = ^ i n several indeterminates X\,.... X n . Namely, 

the orders of power series coefficients equal ord t (Aij ) = a h j. ordt(A) = a and 
the support of X t is Si. More precisely, ([T]) expresses that at least two lowest 
terms of the expansion in power series of the differential equation have the 
same exponents, which is similar to that the tropical equations concern the 
lowest terms of the expansions in Puiseux series of algebraic equations. 

We study solvability of a system of tropical linear differential equations 

min{ap? + ai}, 1 < l < k (2) 

where l<*<n,0<j<?’ and for all finite coefficients ap\ a; 6 Z we have 

0 < off}, ai < M. Thus, the bit-size of © is bounded by knr \og 2 (M + 2). 

We say that a solution T\ ,..., T n of ([2]) is minimal if the inequality 
V g^Ti,...,t 71 < Val Sl ,...,s n holds pointwise for any solution Sj,..., S n of (J2J) . 

Note that ([2]) extends tropical linear systems when for all the occurring 
derivatives x^p we have j = 0. Thus, the complexity bound of testing solv¬ 
ability of 02]) in the next theorem generalizes the similar complexity bound of 
solvability of tropical linear systems from m. a. a- 

Theorem 0.1 If a system (dj) of tropical linear differential equations is 
solvable then it has the (unique) minimal solution. There is an algorithm which 
tests solvability of 1) and in case of solvability yields its minimal solution 
within the complexity polynomial in knrM. 
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Note that S C T C Z + iff the inequality Vais > Vair holds pointwise. For 
Si,..S n ,Ti,... ,T n C Z + we have for the pointwise minimum 

Val(s 1 ,...,s n )v(T 1 ,...,T n ) := f / a/ 5lU r 1 ,...s n uT„ = min {Val Sl ,...,s n , Val Tl ,...,T n }■ 

Assume now that (J2]) has a solution Si,...,S n . If s £ S t such that 
s > r then one can replace S t by adding to it all the integers greater than 
s (while keeping Si,.... S n still to be a solution of (J2J )• Therefore, one can 
suppose w.l.o.g. that for every 1 < i < n either Si is finite and moreover 
Si C {0,..., r — 1} or the complement Z + \ S t is hnite. Thus, if we define 
V := Valy(Si,...,s„) where V ranges over all the solutions Si,..., S n of (|2]) then 
\J can be taken over a hnite number of solutions, hence it can be reduced 
to a single solution Ti,..., T n , thereby V = Valr- l .....T n and T { ,..., T n is the 
minimal solution of (J2J) which proves the first statement of Theorem 10.11 

Also we design a polynomial complexity algorithm for solving systems of 
the type (JU) in case of one variable (n = 1). 

Theorem 0.2 There is an algorithm which tests solvability of a system flf) 
of tropical linear differential equations in one variable x and yields its minimal 
solution in case of solvability within the polynomial complexity. More precisely, 
the complexity is bounded by 0{kr\og{rM)). 


1 Bound on the minimal solution of a tropical 
linear differential equation 


Our next goal is to bound the (hnite) complements Z \ 5). For each 1 < i < n 
with a hnite Z \ Si denote by rrii £ Si the minimal element of S t such that 
mt > r. If for some 1 < l < k the inequality min i,j{a^i +Valsi(j), a{\ > M + r 
holds then for every 1 < io < n, 0 < jo < r for which this minimum is attained: 
a i 0 °,i + Val Sio (j 0 ) = min+ Val Si (j), a i} we have Val s io (jo) = m io - jo- 

Consider a graph G which for each hnite Si, 1 < i < n contains a vertex wy 
and for each Si with a hnite complement Z + \ S l contains two vertices Wi, w°°. 
A derivative corresponds to a vertex w°° iff Vals, (j ) = m % — j (provided 
that Z + \ Si is hnite), else x'p corresponds to w*. Also G contains a vertex wq 
to which corresponds every free term ai, 1 < l < k. 

If there are 1 < l < k, 1 < io, ii < n, 0 < jo,ji < r such that 


a i 0 °J + Val Sio (jo) = + Val Sll (ji) 




minfo^ + Val s .(j), a L } < oo 

hi 


(3) 


then we connect in G by an edge vertices which correspond to the derivatives 
and xf^\ Instead of a[ J Q °J + Vals lQ (jo) could be ai, then we consider the 
vertex wq. 
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If a connected component of G contains only vertices of the form w°° with 
rrij, > r then for each wf° from this component we replace m; by rrii — 1, or 
in other terms, augment Sj by m; - 1, preserving so modified Si,.... S n (for 
which we keep the same notation) to be still a solution of ([2]). After that G can 
be modified (we use the same notation for the modified G), and we continue 
this process. Eventually, we arrive to a solution Si,... ,S n whose graph G has 
no connected component satisfying the described property. 

Therefore, each connected component of G contains a vertex of the form 
Wi 0 (or perhaps, w°f with rrii 0 = r ) fulfilling (J3]) for suitable j 0 , l. Then 
Val Si (jo) A r (there is a possibility that instead of a^ Q °j + Vals lQ (jo) we 
consider the free term ai ), hence Vals 1 (j\) < M + r follows from (J3J). Thus, 
for every vertex wj° from this connected component there is a path in G 
of a length at most n — 1 connecting it with a vertex of the form «y 2 (or 
perhaps, wff with = r). Therefore, there is 0 < j < r such that 
rrii — j = Vals^j) < (n — 1 )(M + r) (one can show the latter inequality fol¬ 
lowing along the path and applying the above argument). Thus, we conclude 
with the following lemma. 

Lemma 1.1 For each 1 < i < n for which Z + \ Si is finite the bound 
mi < N (n — 1 )(M + r) + r holds. 

This lemma extends Lemmas 1.2, 2.2 [4] established for tropical linear 
systems. 

2 Algorithm testing solvability and producing 
the minimal solution of a system of tropical 
linear differential equations 

Now we proceed to design an algorithm which tests whether a system (J2]) is 
solvable and if yes then yields its minimal solution T\,... ,T n . The algorithm 
starts with the setting T\ = • • • = T n = {0,..., N} (see Lemma fLTjl . perhaps, 
being not a solution of ([2]), and then modifies Tf,... ,T n recursively while a 
current Ti,... ,T n is not a solution. If eventually a current Tf,..., T n becomes 
a solution then it is the minimal solution. We show by recursion that for 
any solution S±, ..., S n of (|2]) the pointwise inequality Vals 1 ....,s n A Val Ti,...,T n 
holds for a current T\,... ,T n . 

If a current Tf,..., T n is not a solution of (J2]) then two cases can emerge. 
In the first case there exist 1 < io < n, 0 < < r, 1 < l < k such that 

a finite minimum min+ Vair, (j )} < ai is attained at a single pair 
io, jo- Let ValT io (jo ) = s — jo where s G T io is the minimal element of T io 
such that s > jo. The algorithm modifies Tj 0 discarding s from it. The 
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inequality Vais > V alj still holds for any solution Si,... ,S n of ([2D 
since S) C T, 1 < i < n. Note that if s = rrii 0 < N (see Lemma fl.ID then m lQ 
increases by one. If s = N(— m^ 0 ) then S t0 is hnite due to Lemma 11.11 and the 
algorithm discards from the current T 0 all integers p > N. In the second case 
ai is the unique minimum in min,;j{app + Valrfij), a;}, then system 0 has 
no solution. 

The algorithm terminates when either Ti,..., T n is a solution of (]2]h in this 
case Ti,... ,T n is the minimal solution, or the algorithm detects that (J2D has 
no solution. 

When (J2D is homogeneous, i. e. o; = oo, 1 < l < k, system ([2D has a 
solution with all infinite functions Vals i: 1 < i < n. It can happen that 
the algorithm terminates with all T\,... ,T n being void, that means that the 
infinite solution of © is its unique one. 

To bound the complexity of the algorithm observe that it runs at most 
n(N+l) steps because at each step at least one of the current sets Ti,..., T n C 
(0,..., IV} decreases. The cost of each step is polynomial in knr log M (the 
algorithm for every 1 < i < n stores mi, provided that Z + \ T t is hnite, and 
also stores Ti D {0,..., r — 1}). This completes the proof of Theorem 10.11 

3 Polynomial complexity solving systems of 
tropical linear differential equations in one 
variable 

We design a polynomial complexity algorithm for solving a system of tropical 
linear differential equations ([2D hi one variable x. The algorithm basically fol¬ 
lows the algorithm from Theorem 10.11 with a few modifications. In fact, the 
algorithm designed in this Section is a version of the algorithm from Theo¬ 
rem m the modification consists in that its steps are ordered in a special way 
(observe that at each step of the algorithm from Theorem 10.11 there could be 
several choices of an element to be discarded from the current set T). We use 
the notations from Section [21 

First, if there exists s < r from T such that a hnite minimum min_j{ap + 
Vair(j)} is attained at a unique j 0 for some 1 < l < k and it holds Valr(jo) = 
s — j o then the algorithm discards s from T. This is also a step of the algorithm 
from Theorem 10. 1 1 and we refer to it as a step of the finite type. In other words, 
the algorithm designed in this Section has a preference in discarding elements 
s which are less than r. Denote by So > r the minimal element of T D [r, oo). 

Second, let otherwise s 0 be the only candidate to be discarded from T for 
all the equations from (J2D which are not satisfied by T. Then the algorithm 
from Theorem 10.11 would iust discard sq, while the algorithm under description 
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discards from T possibly more elements at one step. 

For each 1 < l < k consider a unique 0 < jo < r (provided that it does 
exist, i. e. the Z-th equation is not satisfied by T) such that V air (jo) = so — jo 
and ap°) + so — jo = min_j{ap^ + Valr(j)}- Take the maximal pi (or the infinity 
when it is not defined) such that 

a [ i o) + s 0 - jo + Pi < <4 J) + Val T (j) (4) 

for any 0 < j < r for which V air (j) = S\ — j for suitable T 3 Si < r. Observe 
that pi > 1. 

Denote by p the maximum of all such pi. Let p = pi 0 for an appropriate 
1 < lo < k. If p — oo then the algorithm discards from T all the elements 
s > r. Else, if p < oo then the algorithm discards all the elements s from T 
such that So < s < So T p. In other words, the algorithm replaces the minimal 
element So of T ft [r, oo) by so + p, we call this step of the algorithm a jump. 
Clearly, the jump replaces p steps of the algorithm from Theorem 10.11 each 
consisting in discarding just one element from T (so, discarding consecutively 
So, s 0 + 1, • • •, s 0 +p — 1) due to the unique monomial a^ 0 ' 1 + x tjo) at which the 
minimum in ([2]) is attained for the Zo - th equation. 

Observe that after a jump either the (new current) T provides a solution of 
(J2]l or the algorithm can execute a step of the finite type because of the choice 
of p, see (J4j), so discards from T some element s < r. 

As in Section [2] the algorithm terminates when either a current T provides 
a solution of (J2J) (being the minimal solution as it was proved in Section [2]) or 
the algorithm exhausts T (which means that T 0 [0, N] = 0, see Sections [U 
[21) • In the latter case if system ([2]) is homogeneous then it has the (unique) 
infinite solution, otherwise a non-homogeneous system has no solutions (again 
similar to Section [21) . 

To estimate the complexity of the algorithm note that after a jump the 
algorithm executes a step of the finite type, i. e. discards from T an element 
s < r. Therefore, the number of steps of the algorithm does not exceed 2r 
taking into the account that the number of steps of the finite type is less or 
equal than r. To bound the jump pobserve that + s 0 —jo+P = —j\ 

for appropriate j \, si < r (cf. (0}). Since j 0 < r we deduce that s 0 +p < 2r+M, 
hence p < r + M. Thus, one can estimate the complexity by 0(kr log(rm)), 
and we complete the proof of Theorem 10.21 
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4 iVP-hardness of solvability of tropical non¬ 
linear differential equations in one variable 

Now generalizing tropical linear differential equations (see the Introduction) 
we consider systems of tropical non-linear differential equations of the form 

min{ap + E *.“} (5) 

” * (<J)6P 

where the coefficients ap £ Z + and the minimum ranges over a certain (finite) 
family of finite multisets P of pairs We view |P| as the degree of the 

monomial a P + S(ij)eP x i- 

Similar to the case of tropical linear differential equations (see the Intro¬ 
duction), we observe that the solvability of (J5]) is necessary for the solvability in 
power series in t of a non-linear differential equation J2{p} ^P'Wu j)<=p Af > = 0 
where ord t (Ap ) = ap. 

We prove that the problem of solvability (with a set S C Z + similar to trop¬ 
ical linear differential equations, see the Introduction) of a system of equations 
of the form ([5]) is VP-hard already in the case of a single variable x. Mention 
that in [9] NP -completeness of the solvability of tropical non-linear systems 
(in several variables) is established. 

We prove iVP-hardness by means of reducing a 3-SAT boolean formula $ 
in n variables y 0 ,, y n -\ (see e. g. 0) to a system P$ of equations of the 
form ([5]) in a single variable x, preserving the property of solvability. 

The system E< j> contains (linear) equations 

min{V 2 - 7+1 \ 0}, 0 < j < 2n — 1 (6) 

These equations mean that the valuation of each even derivative x^\ 0 < j < 
2n — 1 equals either 0 or 1. Also P$ contains (quadratic) equations 

min{P 2j ' ) + x^ +2n \ 1}, 0 < j < n — 1 (7) 

They mean that either Val(x^ 2 A) = 0, Val(x^ +2n ^) = 1 (which corresponds 
to the value ’’true” of the variable y v 0 < j < n — 1 of $) or Val(x^) = 
1, Val{x^ +2 A) = 0 (which corresponds to the value ’’false” of yj, respectively). 
Finally, for each 3-clause of <h, say of the form ^y 3l V y n V ->yj 3 we add to P$ 
a (linear) equation 

min{x (2il+2n) , x (2j2) , x {2js+2n \ 0} (8) 

Clearly, is equivalent to the solvability of the system obtained by uniting 
(EJ), ([Zj) and (JHJ) for all 3-clauses of <f>. Thus, we have proved 

Proposition 4.1 The problem of solvability of systems of tropical non-linear 
differential equations in a single variable is NP-hard. 
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5 Solvability of systems of tropical non-linear 
differential equations is in NP 

Next we prove that the problem of solvability of systems of k tropical non¬ 
linear differential equations of the form ([2]) of degrees \P\ < d fulfilling the 
bounds: 0 < ap < M, 0 < j < r, 1 < i < n (in an arbitrary number 
n of variables) belongs to NP. First, similar to Lemma 11.11 and using the 
notations from Section [H we show that if a system has a solution (with some 
Si,.... S n C Z_)_) then it possesses a sufficiently small solution. 

Substitute the solution into each equation of the form ([5]), then the valua¬ 
tions of some derivatives can equal m l — j , the valuations of all the other 
derivatives consider as being fixed. We treat the system (after this substi¬ 
tution) as an input of the linear programming problem (expressing that the 
minimum in (J5]) is attained at least at two terms) with respect to the inde- 
terminates m; (for all i for which they are defined), and the fixed valuations 
consider as the coefficients of the input. Therefore, this input possesses a so¬ 
lution with rrii bounded (due to Hadamard’s inequality on determinants) by 
N\ ■= n\ ■ ( M + rd) ■ d n . Note that this bound is worse than the bound on N 
established in Lemma [1 . 1 1 for systems of tropical linear differential equations. 

Since in order to give a solution S\,..., S n it suffices just to specify Si D 
[0,7™], 1 < i < n and m* < N { (for i for which it does exist), we get the 
following 

Proposition 5.1 The problem of solvability of systems of tropical non-linear 
differential equations belongs to NP. 

6 Further research 

Similar to tropical linear systems (cf. [T], [2], [3]) it is an open problem, whether 
one can solve system ([2]) of tropical linear differential equations within the com¬ 
plexity polynomial in knr log M (in other words, within the proper polynomial 
complexity)? In Section [3] a polynomial complexity algorithm is designed for 
testing solvability of systems of tropical linear differential equations in one vari¬ 
able (n = 1). Is there a polynomial complexity algorithm for similar systems 
in, say a constant number n > 2 of variables? 

It is known (see 0. a. ffl) that the problem of solvability of systems of 
tropical linear equations is in the complexity class NP D coNP. Does the 
problem of solvability of systems of the type ([2]) of tropical linear differential 
equations belong to coNPl Proposition 15.11 implies that even a more gen¬ 
eral problem of solvability of systems of the type ([5]) of tropical non-linear 
differential equations lies in NP. 








It is proved in (8] the following coincidence for the closure in the euclidean 
topology: Trop(V(I)) = V(Trop(I )) C M n where / C K[X i,..., X n ] is a poly¬ 
nomial ideal over the field K of Puiseux series and V (I) C K n is the variety 
of /. Does there hold an analogue of this coincidence for differential ideals? 
In other words, is it true that for any differential ideal G in n independent 
variables and a family Si, ..., S n C Z + being a solution of the tropical differ¬ 
ential equation Trop(g) for any g E G, there exists a power series solution of 
G whose tropicalization equals Si,, S n ? 

We say that Si,..., S n is a Laurent solution of (J2J) if for every 1 < i < 
n either Si C Z + is as we considered above or Si = {5} is a singleton for 
some negative integer 0 > b E Z. In the latter case Valb(j) = b — j. This 
corresponds to the order of the j- th derivative of a Laurent polynomial of the 
form ct b (l + 0(t)) for a (complex) coefficient c. If all sets among Si,... ,S n 
are negative singletons then the solvability of (J2J) reduces to the solvability of 
a tropical linear system. The question is, what is the complexity of testing 
whether ([2]) has a Laurent solution? Actually, one can extend this setting 
from Laurent solutions to solutions of the form S t = {b} where b E M \ Z + . 
This corresponds to a necessary condition of solvability of a system of linear 
differential equations in Puiseux series (when b E Q) or in Hahn series (when 
b E M, see e. g. [5]). 

For a tropical linear differential monomial a + x^\ a,j G Z + define its 
derivative as min{o — 1 + x^\ a + x^ +1 ^} when a > 1 or as x^ +1 ^ when 
a = 0 (which mimics the usual derivation law). We spread this definition of 
the derivative to all tropical linear differential equations of the form (J2]) by 
the tropical linearity. The tropical ideal generated by the derivatives of all 
the orders of tropical linear differential equations is called the tropical linear 
differential ideal generated by these equations. Is it possible to test solvability 
of a tropical linear differential ideal? Lest there would be a misunderstanding, 
we note that a solution of a tropical linear differential equation is not necessary 
a solution of the tropical ideal generated by this equation. 

Acknowledgements. The authors is grateful to the Max-Planck Institut 
fiir Mathematik, Bonn for its hospitality during writing this paper. 
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